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Abstract 

The two-dimensional diffusive dynamics of test particles in a random electromagnetic field is 
studied. The synthetic electromagnetic fluctuations are generated through randomly placed mag- 
netised "clouds" oscillating with a frequency oo. We investigate the mean square displacements 
of particles in both position and velocity spaces. As uo increases the particles undergo standard 
(Brownian-like) motion, anomalous diffusion and ballistic motion in position space. Although in 
general the diffusion properties in velocity space are not trivially related to those in position space, 
we find that energization is present only when particles display anomalous diffusion in position 
space. The anomalous character of the diffusion is only in the non-standard values of the scaling 
exponents while the process is Gaussian. 
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I. INTRODUCTION 

Since the paper by Fermi on the origin of cosmic radiation stochastic acceleration 
of charged particles in random electromagnetic fields is a topic of primary interest because 
it plays a crucial role in understanding a lot of dynamical processes in astrophysics, space 
physics and laboratory plasmas (among the huge literature cf. e.g. Ref.s S, li, 3, Q, Q] 
and references therein). In the non-relativistic case stochastic acceleration can be viewed 
as a diffusion process in position and velocity space and it can be studied using a dynamical 



system approach 
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12|. In this framework, the diffusive motion of particles in 



turbulent magnetic fields reveals a rich class of nonlinear phenomena, but the influence of 
electric fluctuations is not well understood, at least not at the same level as for magnetic 
fluctuations. When electric fluctuations are present, the electromagnetic field is not station- 
ary, diffusion in velocity space can occur also in 2-D and it is at least as interesting as the 
usual diffusion in position space. In particular the well known different diffusive regimes in 
position space, namely standard (i.e. brownian-like) and anomalous diffusion, are related 
in a nontrivial way to diffusive regimes in velocity space [lo|, [l^ . In the present paper we 
investigate the diffusive dynamics of test particles subjected to a stochastic electromagnetic 
field. Electric and magnetic fluctuations are obtained through a random distribution of 
oscillating magnetised "clouds" . Different diffusive regimes are found as the oscillation fre- 
quency of the clouds changes, and a nontrivial relationship between diffusive properties in 
position and velocity spaces is observed. 



II. THE STOCHASTIC ELECTROMAGNETIC FIELD MODEL 

We study the non-relativistic behavior of charged particles interacting with a stochastic 
electromagnetic field, the motion of each test particle being described by the usual equations 
of motion 

dr 

^ = "' 

^ = ±[E{r,t)+vxBir,t)] . (2) 
dt m 

A simple synthetic model for the electromagnetic fluctuating quantities can be obtained using 
a vector potential A that lies on the {x, y) plane, namely A = {A^i^r, t), Ay{r, t),0). In the 



gauge where the scalar potential is zero, the electric field is given by -E = {Ex,Ey,0) = 
—dA/dt, while the magnetic field lies in the direction perpendicular to the {x,y) plane, 
namely -B = (0, 0, S) = V x A. With these assumptions the particle dynamics becomes 
bidimensional because the z component of the Lorentz force is zero. 

We consider a model where the electromagnetic fiuctuations are due to the superposition 
of electric and magnetic fields generated by magnetised "clouds" which evolve in time ac- 
cording to some prescribed laws. While in the original model by Fermi cosmic rays interact 
with clouds only through usual collisions , [lljl, here we investigate the case where the 



Lorenz force acts continuously on a test particle. We consider an L x L square as the basic 
cell containing "clouds" . We suppose that the basic cell is repeated in space and that the 
two components of the vector potential are thus constructed through the superposition 
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where ■?/' = e~'^", = |^ — ''"n(i)|/-R, fnif) = {xn{t) , yn{t)) are the coordinates of the n-th 
cloud , R represents the typical spatial extension of the potential produced by a single cloud, 
and = \ f — 'f'mit)\/ R, with being the coordinates of the obstacles in the neighbours 
of the basic cell. According to the hypothesis that the other cells are duplicates of the basic 
cell, we have that = fn + {'>'L,jL) where i and j are integer numbers and at least one of 
them is non zero. Moreover, in order to simplify the sum in Eq. ([3]), R and have been 
chosen in such a way that R is typically smaller than the distance between two nearby clouds, 
therefore in the summation J2m of Eq. (j3]) it is sufficient to consider only the eight nearest 
neighours of the basic cell (i.e. = (±1, 0), (0, ±1), (±1, ±1), (=f1, ±1) )• Since we are 

interested in investigating the diffusive properties of the model and not in reproducing the 
detailed properties of cosmic rays, at variance with the model by Fermi [l| where clouds are 
allowed to move with a given bulk speed, here the motion of the obstacles has been chosen 
to be sinusoidal along both x and y. That is, the coordinates of the n-th obstacle Xnit) and 
yn{t) are given by 

Xn{t) = XnO + an COs[uJ + (4) 
Vnit) = VnO + bn Sm[uJnt + Pn] , (5) 

where Xno and yno are the initial coordinates of the n-th obstacle, a„ and 6„ are the ampli- 
tudes of the oscillations of the n-th obstacle along x and y respectively, Un is the oscillation 



frequency, and a„ and (3n are the initial oscillation phases, randomly chosen within the in- 
terval [0, 2ti]. For the present paper we choose to keep the problem as simple as possible, so 
that we used the same values of motion amplitudes and frequencies for all clouds, namely 
On = bn = a and uJn = ^- The N clouds contained in the basic cell are irregularly distributed, 
that is, with random initial positions (x„o,l/no)- 

The motion equations can be adimensionalised by using the following normalization fac- 
tors: Bq = Aq/L for magnetic field , = {qBo)/m for frequency, to = ujq/^tt for time, 
Vq = Luq for velocity, Eq = AqUq for electric field respectively, so that finally Eq. ([2]) can 
be written as 

— = E{r,t) + vxB{r,t), (6) 

where 

^ Rdin \ |r-r„(t)| J 

I ,x dip \ [x- Xn{t)]dotXn{t) + [y - yn{t)]dotyn{t) \ 

^'<'-''' = ?flae;l r' + ^-^- 

Cj denotes the unit vector along the i-th direction on the plane, and N.T. stands for the 
similar terms given by the clouds in the nearest cells. 

In the stationary limit case {uj = 0) the model is nothing but the electromagnetic refor- 
mulation of the celebrated Lorentz model. The only (marginal) difference is the presence of 
a smooth potential. For uj one has a sort of Lorentz model where the obstacle positions 



are time dependent. In Ref. 
be found. 



^ a dynamical system study in the case of few obstacles can 



III. TEST PARTICLE SIMULATIONS 



The model described above is investigated through test particle numerical simulations. 
The motion equations are solved using a 4th order Runge-Kutta scheme. Particles are 
injected at random positions within the simulation box, with velocities extracted from a 2-D 
Maxwelhan distribution P{vx, Vy) oc exp[— (t>^ + t>^)/2t>^J with f = 3 x 10~^. Since we aim 
to investigate the diffusive properties of the system, all particles have to be tracked for the 
same time interval. The results shown below are obtained by tracking 6 x 10^ particles using 
the parameters = 50, a = 10~^ and R = 0.07, while three different values are used for the 
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FIG. 1: Examples of particle trajectories for u) = 0.01 (left panel) and w = 10 (right panel). 



frequency, namely u = 0.01, oo = 1 and uj = 10. We choose to vary only this last parameter, 
because preliminary results have clearly shown that the model is not significantly sensitive 
to variations of a and R. In Fig. [T]two examples of particles trajectories, for uj = 0.01 and 
u = 10 respectively, are shown. 

In order to study the diffusive properties of the system in both position and velocity 
spaces, the mean square position and velocity displacements are calculated as functions of 
time, that is, {[r{t) ~ r^f) and {[v{t) —vof), where brackets represent averages over the 
particle population. In the case of normal diffusion particles make a Brownian-like motion 
and {[r{t) — r^f) ~ t. On the other hand, if a power law dependence of the mean square 
position displacements on time is found for long times 

{[r{t) - r,f) ~ t^- , (9) 

with Ux 7^ 1/2, the diffusion is said to be anomalous, the cases > 1/2 and < 1/2 being 
called superdiffusion and subdiffusion respectively. Diffusion in the velocity space can be 
studied by investigating the scaling of the velocity mean square displacements with time, 
that is, 

{[v{t) - vof) - t"^^ . (10) 

The mean square displacements of position and velocity are shown in Fig. [2] and Fig. [3] 
respectively for the three vales of u which are used. The best fits with power laws are also 
shown as solid lines. The values of the and Uy exponents obtained from the fits are given 
in Table [H For u = 0.01 the diffusion is nearly normal, since ~ 1/2 and the velocity 
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FIG. 2: Position mean square displacements for u = 0.01 (dotted line), w = 1 (dashed line), 
and a; = 10 (dash-dotted line). The power law best fits are shown as solid lines. The power law 
exponents obtained from these fits are given in Table HI 
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FIG. 3: Velocity mean square displacements for lo = 0.01 (dotted line), w = 1 (dashed line), and 
u; = 10 (dash-dotted line). The power law best fits are shown as solid lines. The power law 
exponents obtained from these fits are given in Table HI 
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TABLE I: Power law exponents Vx ^-nd obtained the from best fits of position and velocity mean 
square displacements, respectively. 





w = 0.01 


UJ = 1 


w = 10 


l^x 


0.49 


0.89 


1.19 




4.7 X 10"^ 


0.42 


0.39 



mean square displacement remains constant after an initial transient. Such an exponent 
value indicates that the system is in a standard diffusion regime. For both uj = 1 and 
= 10 a superdiffusive behavior in position is observed, and the system seems to exhibit 
diffusion also in velocity, with the velocity mean square displacements showing a power 



law dependence. As a consequence o 



exponents Vx and z/„ is not simple 
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the time correlations, the relation between scaling 
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([r(t) - r,f) ^ 2 f r\\v{t')\')C{t\T)dt'dT (11) 



In fact, a straightforward calculation shows that 

ft rt-t' 

lO Jo 

where C(t',T) = {v(t')v(t' + r)) / {\v(t')\^) . In the classical diffusion process, where accelera- 
tion is absent, {\v{t')\'^) is constant, Cit'.r) is independent of t' and {[r{t) — Tq]^) becomes 
proportional to t by integrating over t' . Therefore anomalous diffusion in the position space 
is entirely due to long time tails of C(r). However, in the absence of stationarity the pres- 
ence of C(t',r) makes the scaling exponents Vx and not related in a simple way. In this 



situation only the upper bound z/^; < z/^, + 1 can be derived [lOj]. When the electric field is 
strong, or in other words when very strong time correlations of velocities at different times 
are present, the equality = + 1 (i.e. the maximum allowed difference between and 



u^) seems to occur |10L Il3l |. This indicates that the precise relationship between scaling 
exponents in the position and velocity spaces is strongly dependent on the correlations in- 
troduced in the model. Different models are characterised by different scaling properties. 

n 

For example in the model of Vlahos et al. [7|, where different correlations are used, z/^ ~ 1.15 
is found for position space diffusion, which is similar to the value obtained in the present 
model for uj = 10, but velocity diffusion is weaker, namely z/„ ~ 0.17 with respect ^ 0.39 
(cf. Table H). 

The probability density functions (PDFs) of the two velocity components at different u, 
collected at the end of the numerical simulations, i.e. t = 10^, are shown in Fig. H] together 




FIG. 4: PDFs of the standardised velocity components at u; = 0.01 (dotted lines), to = 1 (dashed 
lines), and w = 10 (dash-dotted lines), collected at the end of the numerical simulations. The 
PDFs at t = are shown as solid lines. 

with the velocity PDFs constructed at t = 0, which are the same in all cases. In order to 
compare the shape of the PDFs for different uj, the PDFs are calculated after transforming 
the variables Vx and Vy into their standardised form by subtracting their averages (vx) and 
(vy) respectively and dividing by the standard deviations cr„^ and a^^ (after this transforma- 
tion the variables have zero mean and unit standard deviation). The shape of the velocity 
PDFs does not change significantly with respect to the initial one, remaining nearly Gaus- 
sian. However, the standard deviations of and Vy increase as u increases (see Table [IT|) . 
i.e. the velocity distributions become broader. 



IV. CONCLUSIONS 



In conclusion we investigated a model of particle diffusion and acceleration in stochas- 
tic electromagnetic fields, focusing on the diffusive properties in both position and velocity 
spaces. As the oscillation frequency of the magnetised "clouds" , which produce the electro- 
magnetic fields, increases, the diffusion in position space goes from a Brownian-like regime 
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TABLE II: Standard deviations of the two velocity components collected at the end of the numerical 
simulations for uj = 0.01, uj = 1, and oj = 10. As a comparison the standard deviations of the 
initial velocities are shown in the second column. 





Initial 


w = 0.01 


tj = 1 


u; = 10 




0.03 


0.03 


0.21 


1.49 




0.03 


0.03 


0.21 


1.49 



to anomalous diffusion, with ballistic motion for the highest investigated frequency. Corre- 
spondingly, diffusion in velocity space also shows a change from Brownian-like behavior to 
anomalous diffusion. However the relation between the diffusion exponents in position and 
velocity does not seem to be simple, depending on the detailed properties of the time corre- 
lation introduced in the model at hand. In the cases where anomalous diffusion is present, 
a broadening of the particle velocity distribution is also found, even if their shape remains 
nearly Gaussian. In other words, the anomalous diffusion mechanism at work is associated 
to an energization of the bulk particle population and not to the existence of high energy 
tails in the velocity distribution. 
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